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Abstract 

We use first-order perturbation theory near the fermionic limit of the (^-function 
Bose gas in one dimension (i.e., a system of weakly interacting fermions) to study 
three situations of physical interest. The calculation is done using a pseudopoten- 
tial which takes the form of a two-body ^'-function interaction. The three cases 
considered are the behavior of the system with a hard wall, with a point where 
the strength of the pseudopotential changes discontinuously, and with a region of 
finite length where the pseudopotential strength is non-zero (this is sometimes used 
as a model for a quantum wire). In all cases, we obtain exact expressions for the 
density to first order in the pseudopotential strength. The asymptotic behaviors 
of the densities are in agreement with the results obtained from bosonization for 
a Tomonaga-Luttinger liquid, namely, an interaction dependent power-law decay 
of the density far from the hard wall, a reflection from the point of discontinuity, 
and transmission resonances for the interacting region of finite length. Our results 
provide a non-trivial verification of the Tomonaga-Luttinger liquid description of 
the 5- function Bose gas near the fermionic limit. 



PACS number: 71.10.Pm 



^^E-mail address: diptiman@cts.iisc.ernet.in 



1 



1. Introduction 

The (^-function Bose gas in one dimension has been studied extensively ever since Lieb 
and Liniger solved it using the Bethe ansatz [1]. However, the wave functions which one 
obtains from the Bethe ansatz are usually quite difficult to work with. It is therefore 
useful to examine other ways of studying the 5-function Bose gas. With this in mind, 
we recently developed a way of doing perturbation theory near the fermionic limit of the 
model [2], which is equivalent to a Fermi gas with weakly attractive interactions. This 
involves the use of a two-body pseudopotential which takes the form of a <5"-function ; 
let us denote the strength of the pseudopotential by a parameter g which will be defined 
below. In Ref. 2, we showed that this perturbative approach correctly reproduces the 
ground state energy up to order g 2 . In this paper, we will apply the pseudopotential 
approach to three situations which cannot be solved using the Bethe ansatz; this will 
illustrate the power of this approach. 

The 5-function Bose gas is an example of a Tomonaga-Luttinger liquid (TLL) [3, 4]. A 
TLL is characterized by two quantities, the velocity of the low-energy excitations v (the 
dispersion relation of these excitations is given by u = v\k\), and the interaction parameter 
K. Once these two parameters are known, the low-energy and long wavelength properties 
of a TLL can be found by the technique of bosonization [3, 4]. In particular, K determines 
the exponents governing the long distance behavior of various correlation functions. We 
will show that the results obtained using the pseudopotential approach are in complete 
agreement with those expected from bosonization. Our study will therefore provide a 
non-trivial check of the pseudopotential approach and will also confirm the expression for 
the parameter K. 

The plan of this paper is as follows. In section 2, we will introduce the 5-function Bose 
gas and the pseudopotential approach for doing perturbation theory near the fermionic 
limit. In section 3, we will consider the 5-function Bose gas with a hard wall. We will 
obtain an expression for the density which is exact to first order in g. We will then find 
the asymptotic behavior of the density far from the hard wall. The expression of the 
density will involve logarithmic factors, which, in section 6, will be recognized as being 
due to an interaction dependent power-law decay of the density which is characteristic of 
a TLL. In section 4, we will consider a model in which the pseudopotential parameter g 
changes disco ntinuously at one point. We will again compute the density exactly to first 
order in g and show that there are oscillations which can be interpreted as being due to a 
reflection of the particles from that point. The amplitudes of reflection from the two sides 
are found to be equal to ±g. In section 5, we will extend the model of section 4 to the 
case of a region of finite length over which the pseudopotential has a non-zero strength. 
Finally, in section 6, we will discuss the TLL approach to the 5-function Bose gas, and 
will compute the Luttinger parameters K and v to first order in g. We will then use 
bosonization to compute the asymptotic behaviors of the density far from a hard wall and 
from a point of discontinuity in the Luttinger parameters. We will then see that these 
results agree precisely with those obtained in the previous sections. We will conclude in 
section 7 by pointing out some possible directions for future research. 
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2. The 5-function Bose gas 



The ^-function Bose gas is defined by the Hamiltonian 

1 d 2 2c 

h = - — y ^ + — y s r Xi - X j ) , (i) 

z '" i l<i<N ux i " L l<i<j<N 

governing N identical bosons moving in one dimension. (We have set h = 1). The 
interaction parameter c has the dimensions of (length)" 1 . We will assume that c is non- 
negative, otherwise the thermodynamic limit (N — > oo) is not well-defined due to the 
presence of bound states with arbitrarily large negative energies. 

We will consider two different boundary conditions for the calculations presented in 
this paper. In one case, the particles will be considered to be on a circle of circumference 
L with the wave functions satisfying periodic boundary conditions. In the other case, the 
particles will be considered to be in a box of length L with the wave functions vanishing if 
any of the particles coordinates is equal to or L (this is called the hard wall condition). 
Since the particles are identical bosons, the wave functions must be completely symmetric. 
There are N\ possible orderings of the particle coordinates, given by < xp 1 < xp 2 < 
• • • < xp N < L, where (Pi, P 2 , ■ ■ •, Pn) is some permutation of the numbers (1, 2, • • •, N). 
If the wave functions are known for any one ordering, say < X\ < x 2 < • • -x^ < L, they 
are known for all other orderings by symmetry. 

It is clear that the model describes noninteracting bosons for c = 0. For c = 00, the 
wave functions vanish whenever any two particle coordinates coincide. We can then carry 
out the unitary transformation ip P — > (— l) p ip P , where ipp denotes the wave function for 
the ordering P, and (— l) p denotes the sign of the permutation P. Under this transforma- 
tion, the wave function becomes completely antisymmetric, i.e., fermionic. Thus c = 00 
denotes a system of noninteracting fermions. (Note that this unitary transformation of 
the wave functions is only allowed if c = 00. At any other value of c, the symmetric wave 
functions do not vanish for x^ = Xj, and the transformation would produce antisymmetric 
wave functions which are discontinuous at those coincident points). Henceforth, we will 
refer to the model with c = 00 as being noninteracting, and the model with c large as 
being weakly interacting in the fermionic sense. 

In Ref. 2, we developed a way of doing perturbation theory near the fermionic limit 
in powers of 1/c. We showed that the perturbation around c = 00 is described by a 
two-particle interaction of the form 

V = - -L £ 5"( Xl - x 3 ) . (2) 

" ic l<i<j<N 

This pseudopotential can be used straightforwardly to first order in 1/c. At second and 
higher orders, some divergences appear which can be cured by a point-splitting prescrip- 
tion [2]. In this paper, we will work only to first order in 1/c and will therefore not 
encounter any divergences. 

We will be interested in the thermodynamic limit in which N, L — > 00 with the density 

N 

Po = T (3) 
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held fixed. In the thermodynamic limit, we find, either from the exact solution [1] or from 
the pseudopotential approach [2], that the ground state energy per particle is given by 

N 6m 3mc 
to first order in 1/c. The chemical potential \x = (dE /dN) L is given by 

AS _ 4^V| 
2m 3mc 

For convenience, we define the Fermi momentum kp and Fermi velocity vp as 

hp = rrp , 

and vp = — . (6) 
m 



We also define a dimensionless parameter 



Equations (4-5) imply that for large values of c, we have a system of fermions with a weak 
attractive interaction of strength —g. 

The energy eigenstates of the system of noninteracting fermions (g = 0) have normal- 
ized wave functions which are given by l/\/iVT times the Slater determinant of a matrix 
M. The entries of M are given by 

M np = ^ n (x p ) , (8) 

where tp n (x) denotes the normalized one-particle wave functions. For the model on a 
circle, ip n is given by 

i> n (x) = -\= e l2 ™ x/L (9) 

where n — 0, ±1, ±2, • • • (if N is odd), and the corresponding one-particle energies are 
E n = (27m) 2 / (2mL 2 ) . For the system in a box with hard walls, the one-particle wave 
functions are given by 

Mx) = sm(— ) (10) 

where n — 1, 2, 3, • • •, and the corresponding energies are _E„ = (im) 2 / {2mL 2 ) . 

Using the pseudopotential given in (2), we can do perturbation theory to first order 
in 1/c as follows. For the noninteracting and weakly interacting systems, let us denote 
the iV-particle wave functions of the energy eigenstates by tyffl (which are normalized to 
unity) and ty n respectively, and the corresponding energies by E^ and E n . To first order 
in 1/c, we have 

* = *«» + y *(°) (jg^jg) (11) 



p(o) p (0) 
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provided that there is no degeneracy at the energy . (This will be true in our calcula- 
tions since we will only apply (11) to the ground state which is unique both on the circle 
and in the box). Note that the norm of fy n in (11) differs from unity only by terms of 
order 1 /c 2 which we are going to ignore. 

The one-particle density p(x) is obtained from the many-particle wave function ^(xi) 

as 

p{x) = N J dx 2 dx 3 ■ ■ ■ dx N , (12) 

where the factor of N has been introduced on the right hand side to ensure that we obtain 
Jq dxp(x) = N. For the system on a circle, we find that the density in the ground state 
of N noninteracting fermions is simply equal to a constant, p(x) = po- For the system in 
a box, the density in the ground state of iV noninteracting fermions is given by 

p(x) = J2j <** 2 (^) • (13) 

n=l L L 

In the thermodynamic limit, we then find that 

. , sin(2/c F x) 

The hole in the density integrates to 1/4, namely, 

f°° 1 

J dx [ p(x) - p (x) ] = - - . (15) 

The oscillations in the density in (14) are caused by reflection from the hard wall at x — 0. 
We will study below how these oscillations are modified by the interactions between the 
particles. 



3. The 5-function Bose gas with a hard wall 

We will now study the behavior of the ^-function Bose gas with a hard wall. We will 
first consider the system placed in a box extending from x = to x = L with hard walls 
at both ends. We will then take the limit L — > oo. Far away from the point x = 0, we 
will see that the leading order oscillatory term in the density has a logarithmic prefactor 
the interpretation of this will be provided in section 6. 

In the ground state of the noninteracting system, the iV particles occupy the 
lowest energy states described by (10) with n — 1, 2, • • •, N. The perturbation in (2) will 
connect this state to two types of states: 

(A) a state of the type ^°\n';n) which differs from the ground state in that only 
one particle is excited from the level n to the level n 1 ', where 1 < n < N and n' > N. 

(B) a state of the type ty(°\n[, n' 2 ; rii, n 2 ) which differs from the ground state ^/q -* in that 
only two particles are excited from the levels ri\ and n 2 to the levels n! x and n' 2 , where 
1 < n 2 < rii < N and n[ > n' 2 > N. 
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The matrix elements of (2) between and states of type A are given by 
V n , n = (¥°\n';n)\V\¥V) = A n , n , 

y-^J _|_ ^ 

where A n i jn = 3(n 2 + ri 2 ) + lOnri if ri + n is even and 1 < — - — < N , 

= - 3(n 2 + n' 2 ) + lOnri if ri + n is even and 1 < < N . 

= otherwise . (16) 

The matrix elements of (2) between and states of type B will be denoted by 
Kii,ji£;ni,ra 2 ; however, we do not need the expressions for these for reasons which will soon 
become clear. 

Following (11), we see that the ground state wave function up to order 1/c is given by 

* = ^ 0) + E E Vn ' n 



n>>N Kn<N 



ri 2 — n' 2 )n 2 1 '(2mL 2 ) 



+ E E * (0) KX 2 ;m,^ V "»un 2 



n >>% N i<n 2 <m<N V " " " " ("? + nl - n? - ri 2 W/(2mL 2 ) " 

(17) 

We can now compute the density using (12). To order 1/c, we see that we get a contri- 
bution from (\I/q ' ) |^>q ' ) ) which is given in (14), and a contribution from the cross-terms 
(^f^\^°)(n'- ri)) arising from the type A states in (17). There is no contribution from the 
cross-terms arising from the type B states since the integration over dx<i ■ ■ ■ dx^ ensures 
that (^/^\^/^(ri l ,ri 2 ;ni,n2)) does not contribute to the density. The vanishing of the 
contribution of type B states leads to a major simplification in the calculations. 

We now find that 

sin(2k F x) 1 ^ . ixrix . nnx A n >. n 

P{X) " P ° + —2^x— = ? Sm( ^ } ri^ ■ (18) 



n' .n 



where A n i- n is given in (16). We now go to the thermodynamic limit. Let us introduce 
the variables k = rcn/L and k' = rrn'/L, so that < k < k F and k' > k F . We replace the 
sums over n and ri by (L/n) J dk and (L/n) J dk' respectively. The condition for non- 
vanishing matrix elements in (16), namely, that ri + n must be an even integer, implies 
that we should put a factor of 1/2 in front of the integrals over k and k' . Further, the 
conditions that (ri ± n)/2 must be an integer lying in the range 1 to N turn into the 
conditions < k! ± k < 2k F . After some manipulations, we find that (18) can be written 
in the form 

, . sm(2k F x) 

1 f kp ji [ 2kF ~ k „ , s\n(kx) s\n(k'x) 



2n 2 c 



rK F rzK F ^K sin (kx sin (kx) r , 9s , , , x 

/ dk / dk' — K J \ — - [ 3(k' 2 + k ) + 10k' k ] . (19) 
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We now define two new variables u — (k' + k)/(2kp) and v — (k' — k)/{2kp). In 
terms of these variables, the limits of the integrals over k and k' in (19) are equiva- 
lent either to Jq du fl^dv or to J 2 dv J^ v _^du. Finally, we can write sin (for) sin(fox) = 
(l/2)[cos(2kFXv) — cos(2kpxu)]. We can now do one of the integrals, either over u or over 
v, to obtain 

, . sin(2foi7'a;) 

— ~ Po9 [ / dv cos(2fop:rt>) {2 — v + — In \v — l|j 
jo ' 2 

fl 1 + u 

+ / du cos(2k F xu) {1 - 2u ln(- )} ] . (20) 

Jo 1 — u 

As shown in the appendix, the total hole in the density is found to be 

1°° I**)-*] = -j - I- (21) 

From (20), we can obtain the asymptotic form of the density at large values of x (i.e., 
for \kpx\ » 1) as shown in the appendix. The density turns out to have an expansion in 
powers of 1/x multiplied by sines, cosines, and logarithms. The leading order terms are 
given by 

/ \ sin(2foirx) r , . r mi \ i cos(2A;ir:r) . . 

p(x) = Po ^^r 1 [1-25 ln^Tre^W) ] + 3g 4^ ' ^ 

where C = 0.5772 ... is Euler's constant. We will see in section 6 that the logarithm is 
a sign that the power of x in the denominator of sin(2fop;r) should really be 1 + 2g (plus 
terms of higher order in g), rather than 1. Namely, 

' L [1 - 2g \n(p x)} (23) 



(p x) 1+2 9 p X 

plus terms of higher order in g. 

4. The 5-function Bose gas with a discontinuity in the interaction strength 

In this section, we will consider another application of the pseudopotential approach. 
We will consider a 5-function Bose gas on a circle, with the interaction parameter in (1) 
being equal to a large value cforO<x<L/2 and 00 for L/2 < x < L. It is clear from 
(5) that the chemical potential will then be different in the two halves of the system; this 
would imply that a state which has the same density everywhere can no longer be close to 
the ground state. We will compensate for this imbalance by adding a one-body potential 
in the region < x < L/2 which is equal to 

5V = . (24) 

3mc 

This ensures that the chemical potential and the density are the same in the two halves 
of the system; therefore the ground states of the noninteracting (c = 00) and interacting 
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(large c) systems are smoothly connected to each other. We will therefore work with the 
following perturbation to the noninteracting system, 

v = ~ -L £ 6 "( x * - ) /(*<) + E /(*0 > 

" ic l<«<j<JV l<i<N 

where /(x) = 1 for 0<i< L/2 , 

= for L/2 < x < L . (25) 

We will then find that there is a reflection from the point of discontinuity at x — 0. 

We will now compute the density of the system to first order in 1/c. In the ground 
state of the noninteracting system, the particles occupy the levels given in(9) with 
n — 0, ±1, ±2, •• •, ±(JV — 1)/2 (assuming that N is odd). As in section 3, the perturbation 
in (25) connects this state to states of types A and B. Type A states differ from 
in that only one particle is excited from a level n lying in the range [— (N — 1)/2, (N — l) /2] 
to a level n' lying outside that range, while type B states differ from ^/q -* in that only two 
particle are excited from levels ri\ and n<i lying in the range [—(N — l)/2, (N — l)/2] to 
levels n[ and n' 2 lying outside that range. After using (11) to write the perturbed ground 
state wave function to order 1/c, we again find that the type B states do not contribute 
to the density. 

The matrix element of (25) between and a type A state ^^(n';n) is given by 

v n , n = (¥°\n';n)\v\¥ 0) ) 

^k F . , n 2 — 1 p 2 L 2 . . £ , 

= T^T-, r [n n + — — — if n - n is odd , 

mcL 2 (n' — n) 12 3 

= otherwise . (26) 

We now go to the thermodynamic limit, and introduce the variables k' = 27in'/L and 
k = 2nn/L. Thus k lies in the range [—kp^kp] while k! lies outside it. We replace the 
sums over ml and n by the integrals (L/2n) J dk' and (L/2n) J dk respectively, and put a 
factor of 1/2 in front of the integrals to take care of the restriction in (26) that n! — n 
must be an odd integer. We then find that the density is given by 

p(x) = po + 4 I"" dk [ dk ' ^T (( f 2 vp M M [k'k-k 2 F ]. (27) 
n 2 J-k F J\k'\>k F (k 12 - k 2 ){k' - k) 

From this expression, it is clear that the choice of the one-body potential in (24) makes 
the matrix element for forward scattering (k' = k = kp or k! = k = —kp) vanish. This 
will lead to the result below that p{x) — po vanishes for x — > ±oo; this is the reason for 
choosing the value of SV as given in (24). 

We now introduce the variables u — (k' + k)/ (2k p) and v — (k' — k)/ (2k p). For v > 0, 
u lies in the range \v — 1| to v + 1, while for v < 0, u lies in the range v — 1 to —\v + 1|. 
After integrating over u, we obtain 

Q f°° , sm(2k F xv) , ^ . „. , . v + 1 . . 

P(x) = Po + f- / dv 1 f ' 2v ■ (v 2 + 1) In | ] . (28) 

2ir Jo v 2 \v — 1 
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From (28), we can obtain the asymptotic behavior of the density as x — > ±00 using the 
methods given in the appendix. We find that 

t v sin(2A; F a;) 
P\x) = Po ~ 9 — ~ for x -> 00 , 

sin(2/ci7:r) „ 

= Po + 9 o for a; -> -00 . 29 

We will now see that the oscillatory terms in (29) can be interpreted as being due 
to reflection of the particles from the point x — 0. Consider a system in which fermions 
coming in from x = 00 get reflected back from x = with an amplitude r+. The density 
at large positive values of x can be shown to be of the form [5, 6] 

p(x) = p + [ r* e~ i2kFX - r + e i2kFX ] . (30) 
On comparing this with (29), we see that 

r+ = -9- (31) 

Similarly, if fermions coming in from x = — 00 get reflected back at x = with an 
amplitude r_, the density at large negative values of x can be shown to be 

p{x) = po + [ r_ e~ i2kFX - r* e i2kpX ] . (32) 
On comparing this with (29), we see that 

r- = 9 ■ (33) 

In section 6, we will see that the values of r + and r_ obtained here agree exactly with 
those obtained from bosonization. 



5. The 5-function Bose gas with a non-zero pseudopotential in a finite region 

We will now extend the analysis of the previous section to the case where the inter- 
action parameter c is non-zero over a finite length /, where / will be held fixed as we go 
to the thermodynamic limit L — > 00. Now there is no particular reason to choose the 
strength of the one-body potential 5V in the finite region < x < I to be the same as 
in (24); as long as 5V is small, the ground state of this system is smoothly connected to 
that of the noninteracting system in which 1/c and SV are zero everywhere. At the end 
of our analysis, we will discover that the reflection from the finite region is zero, i.e., there 
is a resonance in the transmission, if either / has certain special values (given by Tin/ hp), 
or if 5V has a particular value. 

We will work with the following perturbation to the noninteracting (c = 00) system, 
V = - — £ 6"( x t - x 3 ) f(xi) + 5V £ f(xi) , 

1710 ±<i<j<N l<i<N 

where f(x) = 1 for < x < c , 

= for c < x < L . (34) 
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As in section 4, we will compute the density of the system to first order in 1/c. Once again, 
we find that the perturbation in (34) connects the ground state of the noninteracting 
system to states of types A and B, and the type B states do not contribute to the density 
to order 1/c. 

The matrix element of (34) between and a type A state ^°\n';n) is given by 



(¥°\n';n)\V\^ 



(oh 



i4:k F , 1 , , N 2 — I, L 2 5V , r r .„ , ,1 



[ — (n'n + ) - — — ] [ exp{i2vr(n - n') — } - 1 ] . 



L 2 (n'-n) mc 12 47r 2 p 

(35) 

We introduce the variables fc = 2ixn/ L lying in the range [— fop, fop] and and k' = 2nn'/L 
lying outside that range. After going to the thermodynamic limit, we find that the density 
is given by 

k F f , sin{(/c' - k)(x - I)} - sm{(k' - k)x} 



p(x) = po - / dk [ dk' [ 

7l Z J~k F J\k'\>k F 



k F J\k'\>k F (k' 2 - k 2 )(k' - k) 

x[k'k + ^-^]. (36) 
3 Po 

After introducing the variables u = (k' + k)/{2k F ) and v = (k' — k)/(2k F ), and integrating 
over u, we obtain 



g roo 

p{x) = po - — J dv [ 



sin{(2k F (x — l)v} — sin{2k F xv } 



v 2 



r ^ / 9 1 mc5V . , . v + 1 . , 

x 1 2v - <" - 3 + 1,1 W 1 ' (37) 

Using the methods given in the appendix, we find that the asymptotic behavior of the 
density as x — > ±oo is given by 



,2p m5V . sin{2A;p(x — /)} — sin{2/cpx} 



2p m5V sin{2/cp(:r - /)} - sin{2/cp:r} 
= P °- ( ^ + ^pI ) 4^ for s - -oo . (38) 

As in section 4, the oscillatory terms in (38) can be interpreted as being due to 
reflection of the particles from the region < x < I. On comparing (38) with the 
expressions in (30) and (32), we find that 



r + = - i ( ^ + ^ ) sin(M) e"*' 1 , (39) 

OC 7T Po 



and 

2po + N A ^ 

3c 7T 2 Po 



r_ = - i ( ^ + ) sin(M) e iM . (40) 
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We observe that the reflection amplitudes are zero (and the transmission probabilities are 
therefore equal to 1) if either 

where n is an integer, or 

SV = - ^-^ . (42) 
3mc 

We see that up to order 1/c, the resonance condition in (41) is the same as for noninter- 
acting fermions (c = oo); equation (41) is familiar from scattering theory in one-particle 
quantum mechanics. The condition in (42) appears to be new; note that if (42) is satisfied, 
then r± is zero for all values of the length I. From (36), we see that the resonance condi- 
tions in (41-42) imply that the matrix element for backward scattering (k f = —k = kp or 
k' = —k = — /c_f) is zero. 

6. The Tomonaga-Luttinger liquid approach to the 5-function Bose gas 

We will now discuss how the results obtained in sections 3 and 4 may be understood 
using the TLL description of the long wavelength properties of a 5-function Bose gas. As 
mentioned in section 1, a TLL is described by two parameters v and K. These parameters 
can be deduced if one knows the energy of the low- lying states [4]. Let us assume that 
the model is defined on a circle of length L with periodic boundary conditions. Let the 
energy of the ground state be E (N, L). Then one has the relation 

= ^ . (43) 

Next, for a model which is invariant under Galilean transformations (as is true for the 
5-function Bose gas), it turns out that 

vK = v F . (44) 

Since the ground state energy of the 5-function Bose gas is exactly known [1], one can use 
(43-44) to find the values of v and K for any value of the interaction parameter c. From 
(4), we find that 

K = 1 + 2*7, 

v = Dp (1 2g) , (45) 

up to order g. 

The low-energy and long wavelength properties of a one-dimensional fermionic system 
are determined by the modes near the two Fermi points with momenta ±kp- The second- 
quantized fermion fields \1/ can be written in terms of the the fields near the two Fermi 
points as 

■V( x ) = y R ( x ) e ikFX + ^ L (x) e~ ikFX , (46) 

where the subscripts R and L denote right-moving and left-moving respectively. Let us 
define the density operators for these two fields as p R = ^ r ^/r and p F = l. Consider 



11 



a contact interaction of the form [3, 4] 

Vint = J dx[ 92P R (x)p L (x) + y {p 2 R (x) + p 2 L (x)} ] . (47) 
Then one can show that 

v = l{vr + g_2. nvr + B. + 2. )] ». (48) 

Comparing these expressions with (45), we see that for the 5-function Bose gas near the 
fermionic limit, 

^- = 7T- = " 2. • (49) 

We can directly verify the expression for g 2 given in (49). Consider a density-density 
interaction of the form 

Vint = - f f dx dy p{x) U(x-y) p{y) , (50) 



2 

where the density operator follows from (46), 
p{x) = ^(x)^(x) 

= tf R (x)V R (x) + *l(x)* L (x) + ^ R (x)^ L (x)e- i2kFX + ^ L (x)^ R (x)e i2kFX . 

(51) 

If the two-body potential U(x) has a short range, a comparison of (47) and (50) shows 
that the parameter g 2 is related to the Fourier transform of U (x) as 

92 = U(0) - U(2k F ) . (52) 

Comparing (50) and (2), we see that 

U(x) = - — 5"(x) . (53) 
mc 



Then (52) implies that 



92 = L 54 

mc 



which agrees with (49). 

Now we consider the bosonized form of the TLL theory [3, 4]. The bosonic field 4>(x, t) 
is governed by the Lagrangian density 

r - 1 ( d ^ V v ( d(t> \ 2 fnn\ 

L ~ 2vK V dt ) 2ifUJ ■ (55j 

The equations of motion are given by d 2 <p/dt 2 = v 2 d 2 <f)/dx 2 . The excitations of the system 
therefore have the dispersion u = v\k\. 



12 



The technique of bosonization relates the fermi fields ^/r and to exponentials of 
the second- quantized boson field <fi. The exact relationship between the two fields depends 
on the geometrical situation. Let us consider the model of section 3, where the system is 
defined on the half-line x > 0. We denote the incoming (left-moving) fermi field as ^ L 
and the outgoing (right-moving) fermi field as These are not independent fields since 
one is related to the other by reflection at x — 0. One can now "unfold" the half-line to 
the full line and define all the fields to be purely left-moving [3]. The fermi fields $ion 
the full line are related to those on the half-line as 

V L (x) = V L (x) , 
*l(~x) = Vr(x) , (56) 

where x > 0. The bosonized form of this fermionic theory also contains only a left-moving 
boson field 4> L defined on the full line; the two fields are related as 

^ L (x) ~ exp [i {<t> L {x) +<j) L (-x)} + i V^K {<t> L {x) - 4>l{-x)}} . (57) 

Now we can compute the fermion density which is equal to an expectation value in the 
ground state, p(x) = (^(x)^(x)) , where x > 0. Following (51) and (56), this is given by 

p(x) = (& L (x)V L (x)) + (*l(-x)* L (-x)) 

+ (¥ L (-x)^ L (x)} e~ l2kFX + (¥ L (x)^ L (-x)} e l2kFX . (58) 

On using the bosonization expression (57), we find that the first two terms on the right 
hand side of (58) are independent of x; they give rise to a constant which is p . The last 
two terms in (58) give 

(exp [±t2V^K{M-x)-Mx)}}} ~ \- (59) 
Including the factors of e ±l2kFX in (58), we find that 

p(x) - po K -f^ . (60) 

Ob 

We thus see that in the presence of a hard wall, the density of a TLL far from the wall 
has an oscillatory piece whose amplitude decays as l/x K . If K is close to 1 as in (45), we 
see that the amplitude has an expansion in powers of g which is given by (1 — 2g \ax)/x. 
This is exactly what we found in section 3. 

We now turn to the model of section 4, where the Luttinger parameters are given by 
(v, K) for x > and by (vp, 1) for x < 0. We can use the Lagrangian density in (55) 
with these parameters to find the equations of motion for the boson fields. The matching 
conditions at x = turn out to be [7] 

(f)( x = 0-,t) = (j)(x = 0+,t) , 
d(p(x,t) \ v / d(j)(x,t) 



vf 



v / u<pyx, i) \ 



Q x J x=o- K ^ dx ' x=o+ 
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We can now consider what happens when a wave is incident from x = — oo. The equations 
of motion give 



(f>(x,t) = exp [ik(x — vpt)] + r_ exp [—ik(x + vpt)} for x < , 

= t- exp [i — — (x — vt)] for x > . (62) 

The matching conditions in (61) now lead to the following expressions for the reflection 
and transmission amplitudes 



t. 



K - 1 

K + l 
2K 

K + l 



Note that current conservation is satisfied since vp{l —r 2 ^) 
for a wave incident from x = oo, we find that 



(63) 



vt 2 and Vp — vK. Similarly, 



t+ 



1-K 

K + l 
2 

K + l 



(64) 



which satisfies current conservation since v(l — r 2 + ) = vpt\. Upon using (45), we see that 
the reflection amplitudes r± obtained here agree with those obtained in section 4. 

We therefore have the interesting result that a discontinuity in the interaction pa- 
rameters is sufficient to cause scattering, even if there is no other scattering mechanism 
(like an impurity) present in the system. The calculations presented in section 4 can be 
viewed as a microscopic derivation of this interesting phenomenon which had earlier been 
obtained only from bosonization [7]. 

The results in section 5 have implications for the subject of transport through a quan- 
tum wire which is sometimes modeled as a TLL of finite length which is bounded on the 
two sides by Fermi liquid leads [7, 8, 9]. In these models, therefore, the Luttinger pa- 
rameters change discontinuously at the contacts between the quantum wire and its leads. 
For the case of two identical impurities in a TLL, it is known that the transmission reso- 
nances are infinitely sharp at zero temperature [9]. Although the model we have studied 
in section 5 has two points of discontinuity in the interaction parameter, rather than two 
impurities, it is possible that the structure of the transmission resonance will be found to 
change significantly if we go up to higher orders in the interaction parameter 1/c. 



6. Discussion 



We have used the pseudopotential approach to study the behavior of a system of 
fermions with weak attractive interactions. The various situations we have considered 
are not solvable by the Bethe ansatz [1]. This is because the Bethe ansatz only works 
in models which are both invariant under translations and have N commuting operators 
including the Hamiltonian. In such systems, the momenta of the N particles ki, k 2 , • • •, 
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are good quantum numbers, and the wave function can be found exactly as a superposition 
of N\ plane waves. In the absence of translation invariance, there are reflections (at a hard 
wall or at a point of discontinuity in the interaction strength) due to which the particle 
momenta are no longer good quantum numbers. The wave functions are then no longer 
a superposition of a finite number of plane waves, and they cannot be found exactly. 

When the Bethe ansatz fails, the pseudopotential approach seems to be the only way 
to obtain exact results near the fermionic limit, although calculational difficulties may 
restrict its use to low orders in perturbation theory. We have shown how exact expressions 
for the density can be obtained in certain situations. While the agreement between our 
results and those obtained from bosonization for the asymptotic behavior of the density is 
satisfying, we should also emphasize that there are relatively few models with interactions 
in which something can be computed at all distances. 

Our methods can be applied to other problems involving weakly interacting fermions 
in one dimension. For instance, one can study the Kane-Fisher model of a single impurity 
placed in a TLL [9], and the effect of a junction of three or more semi-infinite wires [6]. 
While these problems have been studied earlier using bosonization (and other methods 
which are only valid at long wavelengths [5]), it may be interesting to apply the pseu- 
dopotential method to these situations since we may be able to obtain expressions for 
certain quantities which are valid at all distances. Finally, one can use the pseudopoten- 
tial method to study dynamical quantities like the conductance of a finite length TLL at 
finite frequencies. 

After this work was completed, we found a paper which discusses some properties of the 
one-dimensional Bose-Hubbard model at low densities [10]; the continuum Hamiltonian 
which governs that system is essentially the same as the one studied by us. 
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Appendix 

In this appendix, we will discuss some methods for computing the various expressions 
for the density presented above. 

Let us start with the expression for the hole in (21). Suppose that we have a function 
f(x) defined by the integral 



03(0911)/00/EMR-II. 



f(x) — du cos(axu) h(u) , 



Jo 



(65) 



where h(u) is finite and continuous at u — 0. Then we find 




dx f(x) = lim, 




.2 
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= lim^o I r du h(u) e - 2 « 2 /(4a) ; 

Z V ft JO 



7T 



2a MO) • (66) 

On applying this to (20), we obtain (21). 

Next, we discuss how to obtain asymptotic expressions (for ax » 1) for functions of 
the type 

rug 

fix) = / du cos(axu) h(u) , 
jo 
ruo 

g(x) — du sin(axw) h(u) . (67) 

jo 

If h(u) is finite and continuous for all values of u in the range [0, uq], then the leading order 
expressions for (67) are of order 1/x, and they are obtained by integrating the functions 
cos(axu) and sm(axu). Namely, 

f(x) = — [ h(u ) sin(a:rwo) ] , 
ax 

g(x) = — [ h(0) — h(u ) cos(axu ) } , (68) 

(XOb 

plus terms of order 1/x 2 and higher. These formulae are valid even if m = oo, provided 
that h(oo) = 0. 

Now suppose that the function h(u) has a logarithmic divergence at one point. Some 
examples of this are the integrals [11] 

/ du sin(aa;-u) lnu = — — [ ln(ax) + C ] , 
Jo ax 

1 7T 

du cos(axu) lnu = — , (69) 

2ax 



o 



plus terms of order 1/x 2 and higher. (Here C = 0.5722. . . is Euler's constant). From 
(69), one can show that 

f 1 1 7T 

/ du sin(aa;-u) ln(l — u) = — [ In (ax) cos (ax) + C cos (ax) — — sin(aa;) ] , 
Jo ax 2 

f 1 1 7T 

/ du cos(axu) ln(l — u) = — — [ ln(ax) sin(ax) + C sin(ax) + — cos(ax) ] , 
jo ax 2 

(70) 

plus terms of order 1/x 2 and higher. 

Now consider a function h(u) in (67) which contains a logarithmic divergence at u — 
of the form b lnu, where b is some constant. Then we define a new function h(u) = 
h(u) — b lnu which is finite for all values of u in the range of integration. We can then 
integrate over h(u) using (67-68), and over b lnu using (69). Combining the two gives the 
result of integrating over h(u). Similarly, we can compute the integrals in (67) if h(u) has 
a logarithmic divergence at u — u . Finally, we can also compute integrals over functions 
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h(u) in (67) which have a logarithmic divergence at a point U\ which lies inside the range 
[0, Mo]- We simply divide the integrals in (67) into two parts, one over the range [0,-ui] 
and the other over the range [ui,tio]- Then each of these integrals can be computed as 
explained above. Using all these methods, one can derive the expressions in (22) and (29). 
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